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The Effect of Dust Extinction on the Observed Properties of 
Galaxies in the Near-Infrared. 
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ABSTRACT 

Galaxies behind the Milky Way suffer size reduction and dimming due to their obscuration 
by dust in the disk of our Galaxy. The degree of obscuration is wavelength dependent. It 
decreases towards longer wavelengths. Compared to the optical, the Near InfraRed (NIR) 
Ks band extinction is only « 10% that of the B band. This makes NIR surveys well suited 
for galaxy surveys close to the Galactic Plane where extinction is severe. 

While Galactic obscuration is less prominent in the NIR it is not negligible. In this paper 
we derive empirical relations to correct isophotal radii and magnitudes of galaxies observed 
in the NIR for foreground absorption. We simulate extinction in the J, H and K g bands on 64 
(unobscured) galaxies from the 2MASS Large Galaxy Atlas dJarrett et alj|2003h . We propose 
two methods for the extinction correction, the first is optimized to provide the most accurate 
correction and the second provides a convenient statistical correction that works adequately 
in lower extinction regions. The optimized correction utilizes the galaxy surface brightness, 
either the disk central surface brightness, fig, or the combined disk plus bulge central sur- 
face brightness, elliptical and disk/spiral Hubble types. A detailed comparison between the 
different methods and their accuracy is provided. 

Key words: galaxies: photometry - infrared:galaxies - dust, extinction. 



1 INTRODUCTION 

The effect of dust extinction on point-like objects, stars, is lin- 
early related to extinction. If, for example, the observed mag- 
nitude of a star is Ks = 6™0 and it suffers an extinction of 
Ak^ = 1™5, its intrinsic magnitude will be K° = 4™ 5. Ex- 
tended objects, like galaxies, suffer a further dimming due to the 
loss of their fainter outer regions in the sky background. An ob- 
scured galaxy appears smaller and fainter in the sky than it really 
is. Corrections fo r this ext ra dimming are no n-lin ear as was shown 
bv lFisher&fu iiv (1981); Hauschildt (1987|) and lCameroi] ( 119901) 
in the optical. INagavama et al, ( i2004.) derived corrections to the 
Ks band isophotal magnitudes. Their extinction correction study 
was initiated to correct the isophotal magnitudes for galaxies ob- 
served in the vicinity of the giant radio galaxy PKS 1343-601 cen- 
tred at(i,fe) — (309?7, 1? 7). The average extinction in their region 
was 1'"T From their study they showed that the isophotal 

magnitude correction can be approximated by a linear relation for 
extinction levels Ak, ^ 1™0. Their work also showed that the 
different morphological types of galaxies are affected differently 
by extinction. Galaxies that have an exponential light profile bulge 
(late type) need larger corrections than those with a de Vaucouleurs 
profile (early type), or those with a bulge+disk light profile. 

Statistical results and conclusions derived from magnitude- 
limited and radius-limited galaxy catalogues are unreliable if they 
are not corrected for these extinction effects. Applying extinction 
corrections to galaxy catalogues like the 2MASS Extended Source 



Catalog (2MASX) djarrett et alj200ol) will imply fainter complete- 
ness magnitudes than the currently quoted values. Applying extinc- 
tion corrections to optical catalogues for galaxies observed in the 
Zone of Avoidance (ZOA) helped reduce the optical ZOA from the 
region with extinction < l'." 0, to the region with As < 3'"0 
dKraan-Korteweg & Lahavl2000l) . This reduction was a result of in- 
cluding those galaxies that did not make it into the diameter limited 
catalogues if they were not corrected for extinction. 

Extinction magnitude corrections are crucial when relating 
spiral galaxy magnitudes to their velocity width for the Tully- 
Fisher (TF) relation. The here derived NIR extinction corrections 
will prove inva luable for e.g. the on going whole-sky 2MASS Tully 
Fisher Survey ( [Masters et. alll200^ . especially for galaxies found 
in the ZOA. Application of NIR magnitude extincti on corrections 
will also improve photometric redshift estimations ( |jarrettll2004l) 
for galaxies that are observed in the ZOA. 

Motivated by an ongoing NIR imaging survey to map the 
crossing of the Great Attractor wall across the ZOA where extinc- 
tion levels are severe, we initiated a study of the effect of Galactic 
extinction on galaxies imaged in the J(1.2/im), H{1.6fim) and 
Ks{2.2fim) bands. This will allow us to correct their observed 
isophotal-radii and magnitudes for extinction. In our current ZOA 
survey along the Norma wall (Riad et al., in preparation) we no- 
ticed, that « 85% of the galaxies in our sample were found in 
regions with an extinction Aj < 1™0, while « 14% were found 
in regions with an extinction 1™0 < Aj < 3™0. Only 1% of the 
galaxies were found in regions with an extinction Aj > 3™0. For 
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the H and Ks band ~ 94% and 98% respectively of the galaxies 
were found in regions with Ah,Ks J; 1™0, while « 5% and 1% 
of the galaxies were obscured by 1™0 < Ah,k, < 3™0. Only 
a handfuU of the galaxies were found in regions with Ah,Ks > 
3'."0. We therefore limited the simulation to the extinction range 
Aj^H,Kg = 0™0 — 3™0. In this work and for the purpose of esti- 
mating the extinction suffered by galaxies in our ZOA survey we 
used the general extinction law, 



E{B - V) 



Where E(B — V) is the colour reddening derived from the 
ISchlegel etaLl ( ll998h reddening maps and Ay is th e extinction in 
the o ptical V band. A typical value for Rv is 3.1 jCardelli et al.l 
Il989l) . Extinction in the NIR J, H and Ks passbands w as derived 
using the parametrization given bv lCardellier^ ( ll989h . see Eqs. 

□-El 



Aj = 0.874S(B- 1/), 
Ah = 0.549£(B-V), 
Ak, = 0.365i;(B- V). 



(1) 
(2) 

(3) 



In this paper, we derive NIR corrections to isophotal magni- 
tudes and radii for galaxies obscured by the Milky Way. Section 
2 describes the data set and method. We describe two methods to 
apply the corrections (Sects. 2.2 and 2.3). In Sect. 3 we provide a 
brief comparison between the different methods, and discuss their 
respective reliability. 



2 DATA AND METHOD 

We sele cted 64 galaxies f rom the 2MASS Large Galaxy Atlas 
(LGA) jjarrett et alj l2003h to simulate the effect of extinction 
on isophotal-radii and magnitude for galaxies observed in the 
NIR. 2MASS is an All-Sky NIR survey in the bands J(1.2 ^im), 
H{1.Q iim) and Ks{2.2 iim). Galaxies were selected in such a 
way that they are minimally affected by contamination from neigh- 
bouring sources, and give a fair representation of all morphologi- 
cal types. The sample includes 25 elliptical and lenticular (E/SO), 
and 39 spiral (S) galaxies. Half of the galaxies in the spiral galaxy 
sample are barred with some of them having ring features. The se- 
lected galaxies cover a wide range in galaxy size and brightness. 
The apparent radii range between r2o = 34" — 1614" where r2o 
is the isophotal radius at the surface brightness level /ia' = 20 
mag/arcsec^ . The apparent isophotal magnitude range covered by 
our sample is 7^20 = 1™04 — 1(J"77, where K20 is the integrated 
magnitude within r2o- The galaxies in our sample suffered mini- 
mal obscuration ranging between Ak^ = 0™01 — 0^07, see last 
column in Table[T] The list of selected galaxies is given in Table[T] 
This table lists the most common name of the galaxy, followed by 
the morphological type, r2o, K20, inclination ib/a) and position 
angle (PA) in the Ks band. All the data are taken from ljarrett et al.l 
( I2OO3I). In the table we also give the disk central surface brightness 
Un lFreemanlll97(]|) and the central surface brightness in the Ks 
band. The last column lists the mean exti nction Ak, ■ E xtinction in 
the Ks band was derived from the Schleg el et alj ( Il998l) reddening 
maps and Eqn. [3] The surface brightness profiles for the galaxies 
are taken from the LGjAQ (given in tabular format). 
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^ The surface brightness profile for galaxies in the LGA are found in 
|http://irsa.ipac.caltech.edu/applications/2MASS/LGA/| 



Figure 1. The effect of extinction on the spiral galaxy M88. Left panel: 
surface brightness profile for M88 in the Ks band, with various levels of 
the simulated extinction in the range = CK"0 to 3™0. Right panel; 

integrated isophotal magnitude rriiso with extinction with the same levels 
of the simulated extinction. 



The method we implement here for simulating the effect of 
extincti on on galaxies cl osely follows the precepts of earlier work 
done bv lCameronl l ll99(]h in the B-band. 

It is difficult to define a reliable total radius and magnitude for 
a galaxy jjarrett et~al .l2003h . but more straightforward to work with 
isophotal radii and magnitudes out to a given radius. For this rea- 
son we only consider here connections to isophotal radii and the flux 
within that respective radius. In this paper, we define the limiting 
isophotal radius R as the radii to an isophote of 21.4 mag/arcsec^ 
in the J, 20.6 mag/arcsec'^ in the H, and 20 mag/arcsec'^ in the Ks 
band. Those value s are approximately the la NIR sky background 
jjarrett et alj2003l) . We also define the integrated isophotal magni- 
tude rriiso as the integrated magnitude out to this radius R. 



2.1 Simulating extinction 

The simulation of the effect of Galactic extinction on the appar- 
ent properties of a galaxy was achieved by the inward displace- 
ment of the limiting isophote along the surface brightness pro- 
file. The limiting isophote was shifted in non-uniform steps along 
the galaxy light profile with each step being equivalent to a cer- 
tain extinction level in the investigated band. As an example, the 
left panel of Fig. [T] shows the effect of extinction on the sur- 
face brightness profile of the spiral galaxy M88 in the Ks band. 
The panel indicates increasing levels of extinction at Ak^ = 
0'."0, 0'."7, 1°'4, 1"7, lY^g, T2, 2?^6 and Ak, = 3"'0 ap- 
plied to the galaxy. The extinction levels are represented as hori- 
zontal lines on the plots. For each of the horizontal lines, the part 
of the profile lying below the line represents the obscured part of 
the galaxy while the part above is what remains visible at a given 
extinction. The intercept of the line with the ordinate gives the ap- 
parent limiting isophote, while the intercept of the light profile on 
the abscissa shows the size of the obscured galaxy that remains vis- 
ible. The projection of the point of intercept of the line representing 
Aks = 0™0 with the light profile on the abscissa gives the intrinsic 
radius _R°. The additional dimming caused by the loss of the faint 
outer region is demonstrated by the integrated isophotal magnitude 
profile rriiso shown in Fig.[T](right panel). 

The reduced radii and magnitudes were calculated for each 
step and compared to the original values. The quantities calculated 
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Table 1. The sample of galaxies used for the study of extinction effects on the NIR. 



No. 


Galaxy 


Morphology 


'■20 


K20 


b/a 


PA 


MO 












["] 


[mag] 




[°] 


[mag/arcsec'^J 


[mag/arcsec'^] 


[mag] 


' 


LrL46y / 


E6 


123.8 


6.502 


0.63 


+67.5 


16.10 


12.88 


0.011 


2 


M86 


S03/E3 


151.4 


6.283 


0.67 


— 55.0 


16.37 


13.21 


0.01 1 




M60 


E2 


146.6 


5.825 


0.81 


— 72.5 


15.83 


13.06 


0.010 


4 


M59 


E5 


109.0 


6.866 


0.65 


— 15.0 


16. 1 1 


12.67 


0.012 






El 


115.0 


6.347 


0.92 


-57.5 


16.12 


12.93 


0.015 




M105 


El 


108.0 


6.362 


0.85 


+67.5 


15.92 


12.65 


0.009 


7 


JNLiL.3f44 


E4 


88.8 


7.445 


0.62 


+62.5 


16.18 


13.05 


0.015 




INLiL. / ZU 


E5 


90.5 


7.396 


0.55 


-40.0 


16.14 


13.63 


0.006 


9 




E2 


93.7 


7.024 


0.82 


— 87.5 


16.10 


13.29 


0.008 


10 


iNLiL. 14U / 


EG 


100.2 


6.855 


0.95 


+60.0 


16.25 


13.62 


0.025 


" 




E3 


113.0 


6.800 


0.74 


+45.0 


16.32 


13.36 


0.008 


12 


JNLiL.44/3 


E5 


93.7 


7.269 


0.54 


— 85.0 


15.96 


12.89 


0.010 


13 


IN Li L. 4494 


El -2 


91.1 


7.145 


0.87 


— 07.5 


16.15 


13.04 


0.009 


14 


Kjcc A con 


E2 


69.5 


7.915 


0.75 


+92.5 


16.27 


13.67 


0.010 


" 


Mfl^ 1 111 


SO 


35.0 


9.892 


0.56 


-86.5 


16.15 


15.43 


0.010 


16 


JNLiLzilU 


SO 


91.5 


8.565 


0.24 


+47.0 


15.97 


14.99 


0.039 


17 




SO 


42.2 


8.91 1 


0.45 


+37.0 


15.70 


13.46 


0.016 


18 




SO 


61.5 


9.077 


0.24 


— 72.5 


15.75 


14.63 


0.01 1 


19 


JNLiL.i 1 O 


so 


164.6 


5.937 


0.39 


+45.0 


15.74 


12.27 


0.017 


20 


JNLiL4dJd 


E/SO; i 


134.8 


6.628 


0.84 


— 37.5 


16.36 


13.82 


0.010 


21 




E5 


86.0 


7.546 


0.62 


— 17.5 


16.01 


13.58 


0.007 


22 


M49 


tZ/jU(z ) 


179.2 


5.506 


0.81 


-17.5 


16.11 


12.98 


0.008 


23 


M32 


cE2 


147.4 


5. 139 


0.87 


+ 15.0 


14.96 


11.11 


0.023 


24 




r:,+U-l,pec,J>y 


136.0 


5.904 


0.86 


+27.5 


16.05 


13.64 


0.008 


25 






41.6 


10.161 


0.50 


+67.0 


16.71 


15.97 


0.026 


26 


VT/^^^/i O A A* 

JNLiL4Z44 


SA(s)cd 


157.0 


8.1 10 


0.30 


+45.5 


17.49 


16.82 


0.008 


27 


NGC55 


SB(s)m 


273.5 


6.562 


0.30 


— 67.0 


17.12 


17.03 


0.005 


28 


W(~'i" 1 rvTX * 
rs LiL. 1 U / J 


SB(rs)c 


57.7 


9.690 


0.72 


+42.5 


18.70 


16.99 


0.014 


29 


JMj(_24 / 




141.7 


8.180 


0.55 


+01.5 


18.21 


17.01 


0.007 


30 


NGC24 


SA(s)c 


83.0 


9.215 


0.28 


+43.5 


16.97 


17.03 


0.007 


31 


M33 


SA(s)cd 


499.9 


5.38 


0.80 


+50.5 


17.71 


14.10 


0.015 


32 






165.1 


7.686 


0.40 


-15.0 


16.35 


12.83 


0.017 


33 






217.6 


6.587 


0.80 


+ 109.5 


15.76 


12.97 


0.012 


34 


Ml 00 




150.8 


6.810 


0.73 


— 72.5 


17.28 


13.97 


O.OlO 


35 


Ml 06 


C A D/,,\U.T TXTCD 

J)AB(s)D,LlJNbK 


263.9 


5.598 


0.49 


— 20.0 


15.98 


13.39 


0.006 


36 


M51a 




197.5 


5.601 


0.68 


+57.5 


16.38 


13.57 


0.013 


37 


Ml 09 


sAB(rs)DC,LlJNh,K 


142.6 


7.139 


0.56 


+38.0 


17.17 


14.02 


O.OlO 


38 




SA{s)c 


128.9 


7.365 


0.55 


+87.5 


16.41 


14.54 


0.009 


39 


Ml" /ICO 




107.8 


7.133 


0.81 


+05.0 


16.45 


13.81 


0.011 


40 


JNtj(_4oZD 


C A i....\..V..C.1 

aA(r.s)ab,iyz 


214.8 


5.396 


0.57 


— 70.0 


15.24 


12.38 


0.015 


41 


JN(j(_43y4 


iA(sJa,iy i .9 


201.7 


5.009 


0.54 


+87.5 


16.35 


1 1.87 


0.019 


42 


JNtj(_4/4 


SA(s)0 


46.1 


8.815 


0.99 


+45.0 


16.42 


14.05 


0.013 


43 




SB(s)b;pec;.sp 


143.6 


6.861 


0.30 


+35.0 


15.92 


13.84 


0.006 


44 


jNoc /uyu 


SBc?;sp 


136.4 


8.404 


0.26 


— 48.0 


17.13 


13.65 


0.008 


45 


NGC2442 


SAB(.s)bc;pec 


126.4 


7.071 


0.87 


+27.5 


17.04 


13.61 


0.074 


46 




SBO;pec 


124.2 


6.402 


0.99 


+27.5 


15.77 


12.33 


0.013 


47 


JMj(_ I oUo 


CAD/ ..\U^.C.,T 

iAij(s)bc.;iy2 


131.1 


6.732 


0.42 


-37.0 


15.78 


12.77 


0.011 


48 




C A TI/....1I....C,T 

bAB(rs)Dc,by2 


130.8 


6.501 


0.42 


+67.5 


15.25 


12.84 


0.005 


49 


M88 


SA(r.s)b 


155.4 


6.334 


0.44 


-40.0 


15.56 


13.24 


0.014 


50 


NGC4527 


SAB(s)bc 


141.0 


7.021 


0.34 


+69.5 


15.60 


12.89 


0.008 


51 


M63 


SA(rs)bc 


204.2 


5.728 


0.58 


-82.5 


15.76 


12.71 


0.006 


52 


M98 


SA(s)ab,Hll 


172.8 


7.012 


0.34 


+69.5 


16.54 


13.38 


0.013 


53 


NGC47 


SB(rs)bc 


33.8 


9.988 


0.55 


+85.5 


17.12 


15.49 


0.004 


54 


NGC210 


SA(s)b 


54.1 


8.515 


0.52 


-07.5 


15.85 


14.21 


0.008 


55 


NGC628 


SA(s)c 


125.3 


7.187 


0.86 


+87.5 


17.53 


15.47 


0.026 


56 


NGC772 


SA(s)b 


105.2 


7.440 


0.80 


-45.0 


16.64 


14.08 


0.027 


57 


NGC1187 


SB(r)c 


86.1 


8.315 


0.55 


+52.5 


17.18 


14.28 


0.008 


58 


NGC522 


Sc 


68.9 


9.389 


0.14 


+33.0 


15.45 


15.67 


0.032 


59 


NGC891 


SA(s)b 


225.4 


5.994 


0.24 


+25.0 


15.11 


14.55 


0.024 


60 


M31 


SA(s)b 


1614.0 


1.038 


0.56 


+45.0 


14.81 


11.50 


0.023 


61 


NGC5746 


SAB(rs)b?;sp 


162.2 


6.927 


0.26 


-10.0 


15.10 


14.08 


0.015 


62 


M94 


SA(r)ab 


172.3 


5.169 


0.79 


+85.0 


15.50 


11.70 


0.007 


63 


NGCi300 


SB(s)bc 


131.8 


7.896 


0.48 


-79.0 


17.63 


14.26 


0.011 


64 


M65 


SAB(rs)a 


211.6 


6.113 


0.32 


-08.5 


14.95 


13.312 


0.009 



The table lists the moiphological type, r2o , ^^20 , (b/a.) and PA, values are taken from the LGA jjarrett et alj2003l) . 
The table also lists the Ks band fio and fic- Extinction in the Ks band is also listed. Galaxies marked with * are 
those showing the deviated correction trend explained in the next section. 



as a function of extinction were f{R) and Aruiso, defined as. 



nitudes and radii respectively. The different values of f{R) and 
Arriiso corresponding to the simulated extinction values were then 
calculated, and fitted to the empirical relations following the for- 



where m°j,o, rriiso, R° and R are the intrinsic and absorbed mag- 
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Data fit 

Simulated data 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

Af: [mag] 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

All [mag] 



Figure 2. Tlie calculated f{R) (left), and Arriiso (right), with their fitted 
corrections for the spiral galaxy M88 in the Ks band. 



malism given bv lCameronl ( Il99(]h : 



Arriiso 



(4) 

(5) 



where F, v, a, b are the fitting parameters, and A\ the extinction. 
This parameter fitting was done for all galaxies in our sample. It 
was performed separately for all three J, H, Ks bands. 

The resulting corrections show very comparable trends in the 
three bands, due to the similarity of the light profiles of the galaxy 
in the three bands. To demonstrate the effectiveness of our proce- 
dure we therefore reduce the discussion to the results of the Ks 
band only. Figure |2] displays the calculated values of f{R), and 
Arriiso as a function of Ak,, with their fitted relations for the 
galaxy M88 in the Ks band. The ability of Eqns.|4l|5]to fit the sim- 
ulated data varies among the galaxies. Some galaxies show tighter 
fits compared to M88, while others show more dispersion. The vari- 
ation in the quality of the fitting is expected due to the different 
structures imprinted on the surface brightness profile of the galax- 
ies. 





0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Aj^- [mag] 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Af; [mag] 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Aji [mag] 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Ax [mag] 



Figure 3. Simulated conections f{R) and Arriiso- E/SO top panel, S bot- 
tom panel. 



gives tighter estimates to the corrections than using /ic- On the other 
hand, optimizing the corrections using Hc is convenient since it is 
easily measured and is typically found in large catalogues including 
2MASX. A further point to note is that the NIR fic can be used to 
roughly classify the galaxies as early or late-type galaxies jjarrettl 
l20()d) . In Sect. 2.3 we derive a more general correction which can 
be applied when only the galaxies are classified as early or late type 
but neither fio nor /ic is available. 



2.2 Variation among galaxies (optimized corrections) 

In Fig. [3] we show the magnitudes and radii corrections for all the 
E/ SO (top), a nd S galaxies (bo ttom). The curves in the plots clearly 
show that the ICameror] ( Il990l) relations in their present form can 
not fully account for the variation among the galaxies in the two 
respective groups. 

We investigated the origin of these discrepancies by looking 
at trends with e.g. inclination, central surface brightness /ic, and 
the disk central surface brightness /io extrapolated from the fainter 
outer disk. No significant trend was found between the corrections 
and the variation among galaxies inclinations. Most probably a 
much larger sample than the current one is needed to investigate 
the correlation. A clear correlation between fio and fic with the 
corrections was noticed, however. 

It was noted that galaxies with a low central surface brightness 
fic ^ 17.0 mag/arcsec^ in the Ks band, and low surface bright- 
ness galaxies in general require larger corrections. This is expected 
as those galaxies become obscured much more quickly. Galaxies 
that required larger corrections in our sample are M33, NGC24, 
NGC247, NGC1073, NGC55, NGC4244; they are labeled in Fig. 
[3l Even though the corrections were found to depend on /ic, it was 
in fact found that /io shows the better correlation with the correc- 
tions. 

In this section we describe the method to optimize the cor- 
rections using either no or Hc- Optimizing the corrections using no 



2.2.1 Optimized correction based on fio 

It was realized that the derived quantity, /lo, defined as the disk 
central surface brightness extrapolated from the fainter outer re- 
gions of the galaxy correlates well with the deviation among the 
galaxies. The definition of /io correspo nds to the disk c entral sur- 
face brightness for spirals as defined bv lFreemaiJ l ll97(]h . The bet- 
ter correlation of the corrections with fj,o rather than /tc is expected 
since it is always the fainter outer regions of the galaxy that suf- 
fer most from extinction. To accommodate /io into the corrections, 
Eqns.|4]|5]were re- written as. 



Am^soifio) = F(/io) • 



with 



ao exp(/io • ai), 
bo exp(/io • fei), 
Fo exp(/io ■ Fi), 
Vo exp(/io ■ vi). 



(6) 
(7) 



(8) 



Using the parameter products (axb), {F x v) for each galaxy, 
and analyzing how they vary with /to, we derived these parameters 
(flo ,bo, Fo, Vo, ai,bi, Fi, and Vi). The values are listed in Table|2] 
For each galaxy, the value of /to was found from a linear fitting 
to the disk part of the light profile of the galaxy with fix ^ 17.0 
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Table 2. Parameters ao, bo, Fo, Va, ai, 61, Fi, and vi for the ^0 
optimization Eqn.[8] 



Galaxy 



Param. 


J 


H 


Ks 


Qo 


0.0005 


5.2(e-°6) 


0.0002 


ai 


0.3301 


0.6248 


0.3993 


bo 


0.3659 


15.9697 


1.9398 


bi 


0.0704 


-0.1589 


-0.0289 


Fo 


2.8{e-") 


1.5{e-") 


6.7{e-") 


Fi 


1.3090 


1.3912 


1.3309 


Vo 


8.1329 


236.9696 


18.9573 


VI 


-0.0979 


-0.3075 


-0.1561 


ao 


3.3{e-06) 


3.9{e-06) 


0.0001 


01 


0.6093 


0.6203 


0.4406 


bo 


9.7884 


7.0789 


2.6909 


bi 


-0.1116 


-0.095 


-0.0395 


Fo 


8.1{e-") 


1.4{e-i0) 


2.4(e-7) 


Fi 


1.2277 


1.2443 


0.8365 


Vo 


201.5666 


56.9891 


17.5403 


VI 


-0.2693 


-0.2055 


-0.1415 



Table 3. Pai'ameters ao, 60, Fo, Vo, ai, bi, Fi, and vi for the fic 
optimization Egn. llll 



Galaxy 


Param. 


J 


H 


Ks 




do 


1279 


0547 


1598 




ai 


O.UOD J 


0707 


— U.UUU-L 




ft 


8397 


1 3805 


7292 


E/SO 




0248 


—0.0115 


0381 




Fo 


0.0073 


0.0034 


0.0151 




Fi 


0.1889 


0.2574 


0.1632 




Vo 


1.4504 


2.1299 


0.8602 




VI 


0.0046 


-0.0253 


0.0420 




ao 


0.0068 


0.0095 


0.0157 




01 


0.1953 


0.1824 


0.1641 




bo 


2.7721 


2.5213 


1.3653 


S 


bi 


-0.0455 


-0.0397 


0.0009 




Fo 


0.0003 


0.0007 


0.0012 




Fi 


0.4136 


0.3833 


0.3751 




Vo 


6.4910 


4.6215 


2.3296 




VI 


-0.0866 


-0.0674 


-0.0267 



mag/arcsec . The intercept of the straight line fitted to the disk part 
of the light profile with the surface brightness axis was then ^0 ■ The 
disk part of the light profile was identified by visual inspection. 

The goodness of f{R, ^lo) and Arriiso in describing the cor- 
rection curves for the different galaxies accurately was found to 
be very sensitive to the value of ^o- This shows the importance of 
having an accurate light profile of the galaxies, and hence /lo, to be 
able to use the ^0 optimized correction. 

2.2.2 Optimized correction based on /ic 

For optimizing the extinction corrections using the central surface 
brightness /ic we used the same methodology as for the /io opti- 
mization. To incorporate in the corrections we re-wrote Eqns.|4] 
and|5]as, 



with 



a{fj,c) 



ao exTp(nc ■ ai), 
bo exp(/ic ■ bi), 
Fo exp(/ic ■ Fi), 
Vo exp(/j,c ■ vi). 



(9) 
(10) 



(11) 



Again we determine the products (a x b) and {F x v) and 
how they vary with /ic and hence derived the parameters 
(flo ,bo,Fo,Vo,ai,bi,Fi, and vi). The values of these parameters 
are given in Table |3] 

To quantify the performance of the two optimization correc- 
tions, we computed the difference between the corrections as given 
by the simulated data and as given by the optimized corrections 
using both the parametrization /(i?,/io), Amiso{fJ,o) and using 
f{R, fJ-c), AmisoifJ-c)- The comparison was made for all galaxies 
in the Ks band in the extinction range Ak^ = 0™0 ~ 3'."0. A sum- 
mary of the comparisons is listed in Table|4l and displayed in Fig.|4] 
for the /io optimization, and Table[5] Fig.|5]for the /ic optimization. 

The error bars in Figs. [4] |5] and the plots that follow are the 
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Figure 4. Comparison between simulated corrections and mo optimised 
con'ections derived using f{R, /to) and Amiso (R, t^o), top: E/SO, bottom: 
S galaxies. 



standard deviation error defined as SDE — a/VN, with cr, the 
standard deviation of the binned values and N the number of data 
points in the bin, which are 25 for E/SO, and 39 for S. 

The plots and tables show that the accuracy of the corrections 
are extinction dependent, showing smaller deviations and less dis- 
persion at low levels of extinction. It is also evident that spiral 
galaxies show larger deviations and more dispersion than ellipti- 
cal and lenticular galaxies. The relatively larger errors for the spiral 
galaxies can be attributed to the more structural features like spi- 
ral arms, bars traced by their light profiles. The comparison also 
shows that the corrections parametrized using the disk central sur- 
face brightness show less systematic shifts and tighter dispersion 



6 Ihab F. Riad et al. 



Table 4. Comparison between simulated corrections and /iq optimized corrections at tiie extinction levels 
Ak = 0™5, 1™0, 2™0and3™0. 



Galaxy 


Param. 




SDE 


1™0 


SDE 


21" 


SDE 


3™0 


SDE 


E/SO 


f(R) 


0.0012 


0.0071 


0.0017 


0.0130 


0.0087 


0.0236 


0.0424 


0.0433 


E/SO 


Aniiso 


0.0004 


0.0022 


0.0002 


0.0035 


0.0050 


0.0091 


0.0218 


0.0280 


S 


fiR) 


0.0042 


0.0104 


0.0071 


0.0248 


0.0748 


0.1476 


0.7993 


1.1632 


S 


Aniiso 


-0.0016 


0.0069 


-0.0090 


0.0156 


-0.0070 


0.0457 


0.0479 


0.1042 



The ^0 optimized coiTections were calculated using Eqns. |6] [T] The value of for Each galaxy was found by 
performing a lineal' regression to fainter outer region of the galaxy light profile with > 17.0 mag/arcsec^ and 
measuring its intercept with the ordinate. 



Table 5. Comparison between simulated coiTections and optimized conections at the extinction levels 
Ak, = 0™5, 1™0, 2™0and3™0. 



Galaxy 


Param. 


0™5 


SDE 


1™0 


SDE 


2l"0 


SDE 


3™0 


SDE 


E/SO 


f{R) 


0.0010 


0.0080 


0.0025 


0.0170 


0.0148 


0.0437 


0.0623 


0.0972 


E/SO 


Aniiso 


0.0004 


0.0044 


0.0000 


0.0096 


0.0037 


0.0188 


0.0186 


0.0276 


S 


f{R) 


0.0057 


0.0157 


0.0215 


0.0449 


0.2520 


0.2398 


2.0472 


1.4042 


S 


Aniiso 


0.0022 


0.0175 


-0.0011 


0.0388 


0.0252 


0.0868 


0.1356 


0.1451 



The fic optimized corrections were calculated using Eqns.[9l lTol 
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Figure 5. Comparison between simulated corrections and the fic optimised 
corrections derived using f{R, fic) and Aniiso [R, fJ-c), top: E/SO, bottom: 
S galaxies. 



than the corrections using the central surface brightness optimiza- 
tion. The better performace of fic in describing the corrections re- 
sults from the fact that it is the outer galaxy disk that suffers more 
of the obscuration. 

It is also worth noting that « 85%, 94% and 98% of the 
galaxies in the J, H and Ks bands in our ZOA survey are in re- 
gions with an extinction Aj^h,Ks ^ 1™0 where the corrections 
show very little dispersion. 



2.3 Average behaviour (average correction) 

The application of the optimized method to estimate the obscura- 
tion corrections for galaxies requires the knowledge of the value of 
/io and hence the galaxy light profile. In many cases that informa- 
tion is not available, or hard to obtain, like trying to correct isopho- 
tal radii or magnitudes for galaxies in a galaxy catalogue such as 
the 2MASX. For such cases an average correction independent of 
the galaxy light profile is required. In this section we derive such 
an average correction. A comparison between the application of the 
optimized and average corrections is given later in Sect. 3. 

The all elliptical galaxies and the majority of the spiral galax- 
ies (85% of the spirals) show comparable trends (see Fig. |3]l, 
while only a few galaxies (M33, NGC24, NGC247, NGC1073, 
NGC55, NGC4244, marked on the plots) deviate from the aver- 
age behaviour. To produce an average correction for each galaxy 
family E/SO or S, we excluded the strongly deviating galaxies. For 
the ones showing the similar behaviour (25 E/SO and 33 S galax- 
ies) their simulated corrections were binned, resulting in an average 
correction curve for each family in each of the three bands. These 
average curves for the J, H and Ks bands are displayed in Fig. [6] 
For comparison the B-band correctio ns are displayed as well. The 
latter are taken from|c amero ril( ll99(]h . 

As expected, the error bars grow with increasing extinction. It 
should be noted that the error bars are biased with our selection of 
galaxies representing the average behaviour. If we would include 
the strongly deviating galaxies, the error bars would obviously be 
larger. But given that the majority of the galaxies follow the general 
behaviour it would be unreasonable to include those outliers for the 
correction. It is also obvious from Fig.[6]that our correction curves 
are similar to those given by Cameron for the B band, especially 
the shorter J and H bands. The average correction curves were then 
fitted with Eqns.|4][5] The respective fitting parameters are given in 
Tabled 

To gain more insight on the variation of the average correc- 
tion among the galaxies, we calculated the difference between the 
simulated corrections and the average corrections for each galaxy 
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Figure 6. Average correction curves for the 25 E/SO (top), and 33 S galaxies 
(bottom) representing the average correction in each respective galaxy class, 
for the bands J, H, Kg and B band. 



0.20 
0.15 
0.10 
0.05 
0.00 



(e7|o)' 









(eTso) 










14.5 15.0 15.5 16.0 16.5 17.0 



14.5 15.0 15.5 16.0 16.5 17.0 



0. 
0.6 
14 0.4 
5 0.2 
0.0 













I 0.4 

£ 0.2 















.5 15.5 16.5 17.5 18.5 



14.5 15.5 16.5 17.5 18.5 



Figure 7. Difference between simulated corrections and average correction 
for an extinction value of Aj^^ = 1™0. Top: E/SO, radius correction left 
panel, magnitude coiTection right panel. Bottom: S, radius connection left 
panel, magnitude coiTection right panel. The values for the galaxy M88 are 
indicated on the spiral panels with the symbol ®. 



Table 6. Average fitting parameters a, b, F, and v for E/SO and S galaxies. 



Galaxy Param. 



H 





a 





1393 


± 





0001 





1438 


±0 


0001 





1445 


±0 


0003 


E/SO 


b 


1 


2212 







0005 


1 


2056 


±0 


0004 


1 


2555 


±0 


0018 




F 





1111 


± 





0001 





1175 


±0 


0002 





1328 


±0 


0005 




V 


1 


5872 


± 





0008 


1 


5559 


±0 


0019 


1 


5950 


±0 


0041 




a 





1169 


± 





0007 





1185 


±0 


0006 





1453 


±0 


0005 


S 


b 


1 


5650 


± 





0056 


1 


5718 


±0 


0047 


1 


4778 


±0 


0032 




F 





1346 


± 





0008 





1442 


±0 


0007 





1982 


±0 


0009 




V 


1 


9080 


± 





0064 


1 


8755 


±0 


0053 


1 


7093 


±0 


0049 



at Aks ~ 1™0. The average correction was calculated using Eqns. 
mis] and the parameters in Table [6] In Fig.|7]we plot these differ- 
ences against the value of fio for the galaxy in the Ks band. We 
find that E/SO galaxies with fio ^ 16.2 mag/arcsec^ and S galaxies 
with fio ^ 16.6 mag/arcsec^ are generally underestimated by the 
average correction, while brighter galaxies are overestimated by the 
correction. 

To assess the accuracy of the average correction as a function 
of extinction, we made a comparison between the simulated correc- 
tions and those derived from the average correction in the extinc- 
tion range Ak^ = 0™0 - 3™0. The results are plotted in Fig. [8] 
A summary is given in Table [7] where we list the mean difference 
between the simulated data and the average correction value at the 
extinction levels Ak, = C^S, 1™0, X^O and 3'"0, as well as their 
scatter. 

Similar to the optimized corrections, the average correction 
method performs better for lower levels of extinction. Comparing 
Figs. El H] and [S] Tables [4] |5] and [7] we notice that the optimized 
corrections are more accurate than the average correction method. 
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Figure 8. Comparison between simulated corrections and average correc- 
tions. 



3 DISCUSSION 

In this paper we present two methods to correct the isophotal mag- 
nitudes and radii for galaxies observed in the NIR obscured by fore- 
ground extinction. The optimized correction requires knowledge of 
the galaxy's light profile. The use of /i^ to estimate the corrections 
is useful as it can also be independently used to roughly categorize 
the galaxies as early or late type. The average correction method is 
more straightforward. It gives average corrections at each extinc- 
tion level in the J, H, Ks observed wavebands, and only requires 
the classification of a galaxy as early or late type. 
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Table 7. Comparison between simulated corrections and the average corrections at extinction levels 
Ak, = 0™5,1™0, 2™0and3™0. 



Galaxy 


Param. 


0™5 


SDE 


1™0 


SDE 


2™0 


SDE 


3™0 


SDE 


E/SO 


f{R) 


0.0085 


0.0081 


0.0108 


0.0170 


0.0000 


0.0450 


-0.0457 


0.1089 


E/SO 


Aniiso 


0.0425 


0.0046 


0.0045 


0.0108 


-0.0020 


0.0266 


-0.0126 


0.0501 


S 


f{R) 


0.0348 


0.0177 


0.0913 


0.0530 


0.4994 


0.2961 


2.9797 


1.6339 


s 




0.0406 


0.0212 


0.0920 


0.0502 


0.2365 


0.1236 


0.4594 


0.2188 



Average corrections were calculated using Eqns.|4][5] and the parameters in Table|6] 



To compare the accuracy of the average and the /xq method, 
we used their comparison with the simulated corrections, see Figs. 
|4]and[8]and Tables|4]and|7] For spiral galaxies estimating the mag- 
nitude corrections using no optimization, Arrnsoii-io) shows a sys- 
tematic shift A = -0™01 with a SDE = 0™0 for an obscura- 
tion level of Ak^ = 0™5. Meanwhile using the average correction 
shows a A = (r04 with SDE = 0'"02 at the same level of ex- 
tinction. At Aks = 1™0, Amiso(/io) has a systematic shift of 
A = -0™01 and a SDE = 0™02, while the average correction 
has a A = 0'?09 with a SDE = 0™05. Radius corrections for 
spiral galaxies revealed a similar trend for the two methods, (see 
Tables |4]and|7}- The comparisons for the elliptical galaxies are also 
given in Tables|4]and|7] 

Figures [4] [8] and Tables [4] [7] give the comparison 
between the optimized and average corrections as com- 
pared to the simulated corrections for the extinction values 
Ak, = 0™5, 1™0, 2™0, 3'"0. They clearly emphasize that the 
Ho optimized correction is more accurate compared to the average 
correction method. But the average correction remains more useful 
when applying the corrections to galaxy parameters from a galaxy 
catalogue. 

The /ic optimized correction shows larger shifts and more dis- 
persion than the /io optimized correction, but smaller shifts and 
tighter dispersion compared to the average correction. 

In the following we give some average correction values to 
correct magnitudes and radii of obscured galaxies. The average cor- 
rection estimates a 0^13 connection to the isophotal magnitude of 
elliptical galaxies at Ar-^ = 1™0. This magnitude correction is 
over and above the Ak^ ~ 1™0 correction. The magnitude correc- 
tion shows a systematic shift of A = 0™01 with a SDE = O'^^Ol 
at Ak, = 1™0. As a function of radius, ellipticals are estimated 
to be 28.3% smaller in radius at Ak^ ~ 1™0, when using the av- 
erage correction. The radius correction shows a A = 0.6% and 
SDE = 0.9% at the same extinction level (see Table[8}. 

The isophotal magnitudes of spiral galaxies at Ak,, = 1™0 are 
expected to be ff}^20 brighter when applying the average correction. 
They show a systematic shift of A = 0'?^09 with a SDE = (r05 
at Ak, ~ 1™0. The average corrections predict that spiral galax- 
ies appear 28.4% smaller at Ak,, = 1™0. The corrections show a 
systematic shift of A = 4.4% with a SDE = 2.6% at the same 
obscuration level (see Table [8}. The table also lists the expected 
magnitudes and radii corrections respectively at the extinction lev- 
els Aks = 0™5, 1"0 and 2™0. The table also lists the systematic 
shifts and the SDE. The magnitude values listed in Table [8] give 
the additional dimming, the galaxy size reduction, the systematic 
shift when using the average correction and their respective SDE. 
The positive systematic shifts indicate that the average correction 
under estimates the obscuration corrections. 

It is worth mentioning that our corrections agree well with 



corrections given bv lNagavama et all ( |2004|) for galaxies in the Ks 
band obscured by Ak, ^ 1™0. In their work they estimated the ex- 
tra dimming to elliptical galaxies to be Arriiao = 0^15 compared 
to Arriiso = 0™ 13 as expected by our average correction. For spiral 
galaxies they estimated the correction to be Arriiso ~ 0™18 which 
also agree w ell with our expected co rrection of Arriiso = 0™20. 
Compared to iNagavama et al.l ( |2004|) corrections, our results are 
useful for higher extinction levels in the three NIR bands J, H and 



Ks i.e A 



J.H.K, 



3^0. 



4 CONCLUSION 

We present two methods to correct galaxies for extinction in the J, 
H and Ks bands. The optimized correction methods are more ac- 
curate than the average correction. However the average correction 
method is more straightforward to apply as it requires no knowl- 
edge of the light profile of the galaxy but only the classification of 
galaxies as early or late types. The extinction corrections that we 
present here are consi dered as a NIR ex tension to those for the B 
band derived before bv lCameronl jl990l) . 

These corrections will be invaluable to the analysis of large 
scale structures in the ongoing NIR galaxy survey along the Norma 
Wall in the ZOA. It will also be applicable to other galaxy sur- 
veys e.g. 2MASX or prospective ESO galaxy surveys e.g. VISTA 
Kilo-Degree Infrared Galaxy Survey (VIKINcfl) using the VISTA 
telescope. 
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Table 8. Average magnitudes and radii corrections at A^^ = CK"5, 1™0, 2™0 3™0 







0™5 






1™0 






2™0 






3™0 




[mag] 


Arriiso 


A 


SDE 


Aniiso 


A 


SDE 


Arriiso 


A 


SDE 


Aniiso 


A 


SDE 


E/SO 
S 


0.044 
0.061 


0.004 
0.041 


0.005 
0.021 


0.133 
0.198 


0.005 
0.092 


0.011 
0.050 


0.401 
0.648 


-0.002 
0.236 


0.027 
0.124 


0.766 
1.296 


-0.013 
0.459 


0.050 
0.219 


[%] 


fim 


A 


SDE 


f{R) 


A 


SDE 




A 


SDE 


/{«) 


A 


SDE 


E/SO 
S 


13.0 
11.3 


0.7 
2.7 


0.6 
1.4 


28.3 
28.4 


0.6 
4.4 


0.9 
2.6 


54.8 
60.6 


0.0 
6.5 


0.9 
4.1 


73.3 
81.7 


-0.3 
6.5 


0.8 
4.3 



The table lists the extra dimming and the radii reduction as estimated by the average corrections. It also lists the systematic 
shift A for comparing the simulated correction and the average connection with the respective SDE. The positive systematic 
shifts indicate that the average connection under estimate the extinction conection. 
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